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1  Statement  of  the  Problem  Studied 

In  the  recent  history  of  the  study  of  the  ABL,  supercomputers  and  computational  fluid 
dynamics  have  replaced  the  traditional  tools  of  Reynolds-decomposition  [1].  In  particular  the 
development  of  the  successful  direct-numerical  simulation  (DNS)  and  large-eddy  simulation 
(LES)  methods  has  made  the  field  of  turbulence  simulation  particularly  accessible.  Each  of 
those  two  methods  can  be  used  to  simulate  ABL  turbulence.  While  DNS  is  an  exact  solution 
to  the  governing  equations,  it  is  limited  by  current  computational  abilities  to  low  Reynolds 
number  turbulence  studies;  a  “Reynolds  number  similarity”  hypothesis  is  used  in  order  to 
extend  it  to  the  high  Reynolds  number  flows  of  practical  importance  [2].  Unlike  DNS,  LES  is 
not  exact.  It  uses  a  hybird  of  simulation  and  modelling  that  filters  the  governing  equations 
into  “subgrid”  and  “resolvable”  scales,  and  models  the  effects  of  the  subgrid  scales  on  the 
resolvable  scale  quantities  [3]. 

While  traditional  turbulence  methods  are  well  grounded  in  experimental  results,  current 
simulation  techniques  are  not;  an  outstanding  question  in  using  DNS  and  LES  is  how  well 
they  compare  with  experimental  atmospheric  data  [1].  Wind  profilers  [4,  5,  6],  sodars  [7,  8], 
FM-CW  radars  [9,  10]  radio  sondes  [11],  and  in  situ  sensors  [12,  13]  have  all  been  used 
with  some  success  in  measuring  the  structure  of  the  ABL.  Each  of  those  methods,  however, 
only  produces  a  one-dimensional  vertical  profile  of  the  atmosphere;  both  DNS  and  LES 
produce  three  dimensional,  time-varying  fields  of  temperature,  velocity,  and  other  quantities 
of  interest.  Experimental  verification  of  those  models  would  be  optimally  accomplished  by 
measuring  three-dimensional,  time-varying  fields  on  scales  comparable  to  those  used  by  DNS 
and  LES. 

This  report  summarizes  MIRSL’s  activities  under  the  Department  of  Defense  University 
Research  Initiative  (URI)  to  develop  a  radar  remote  sensing  system  that  is  able  to  measure 
atmospheric  turbulence  for  ABL  studies,  including  the  verification  of  LES  results.  This 
radar,  called  the  Turbulent  Eddy  Profiler  (TEP)  is  capable  of  imaging  the  structure  of 
turbulence  throughout  a  conical  volume  that  extends  from  the  ground  to  the  top  of  the 
ABL.  Preliminary  field  measurements  made  at  Duck,  NC  and  Rock  Springs,  PA  indicate 
that  TEP  is  able  to  resolve  the  three-dimensional  and  velocity  fields  at  spatial  and 
temporal  scales  comparable  to  LES  computations. 

The  ability  of  TEP  to  measure  fields  is  an  important  feature  of  the  instrument. 
Volume-imaging  lidars  [14,  15]  can  measure  three-dimensional  velocity  fields  on  time  scales 
on  the  order  of  minutes.  The  lidar  senses  the  backscattered  return  from  aerosols,  however, 
and  thus  cannot  measure  directly.  TEP  is  currently  the  only  instrument  to  provide 
volume-images  of  on  the  time-scales  necessary  for  the  study  of  the  evolution  of  turbulent 
structures.  Recent  results  from  our  colleagues  at  PSU  [16,  17]  suggest  that  LES  could  be 
used  to  predict  acoustic  and  electromagnetic  scattering  properties  by  generating  statistics  of 
“local,”  or  pixel-sized  values.  Experimental  verification  of  those  LES  predictions  was  a 
focus  of  this  program. 
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Figure  1:  TEP  field  deployment  site  plan. 

2  Summary  of  Results 

MIRSL  and  PSU  were  awarded  URIs  in  1992  to  further  understanding  of  the  atmospheric 
boundary  layer.  The  UMass  URI  supported  the  creation  of  a  digital  beamforming,  915  MHz 
radar  that  could  remotely  sense  the  index  of  refraction  structure-function  parameter 
through  the  measurement  of  the  volume  backscatter  [18]: 

T)  «  0.38  Cl 

where  A  =  0.328m  is  the  radar  wavelength.  In  addition  to  obtaining  Cl  fields,  this  coherent 
radar  is  also  capable  of  measuring  fields  of  velocity  vectors  and  spectral  widths. 

Construction  of  the  TEP  radar  system  began  in  1992  and  was  completed  in  1996.  A  site 
plan  for  TEP  is  shown  in  Figure  1.  A  48’  trailer  houses  the  transmitter,  radar  electron- 


Figure  2:  TEP  array  assembly  at  Rock  Springs,  PA.  The  TEP  transmitter  antenna  and 
trailer  are  visible  in  hte  background. 

ics,  and  data  acquisition  system  and  is  used  to  transport  the  entire  system.  A  pyramidal 
horn  transmit  antenna  is  fed  by  a  25  kW  peak-power  transmitter  and  illuminates  a  25  de¬ 
gree  conical  volume  above  the  receiver  array.  The  receiver  array  consists  of  60  microstrip 
patch  antennas,  each  with  a  low-noise  receiver,  digital  accumulator,  and  control  and  storage 
electronics.  Radar  echoes  are  coherently  integrated  over  short  intervals  (<  10  ms)  in  the 
accumulator  and  then  stored  on  disk  or  tape.  Figure  2  shows  a  photo  of  the  array  being 
assembled  at  Penn  State  University’s  meteorological  field  site  at  Rock  Springs,  PA. 

Beamforming  is  performed  off-line.  Software  produces  focused  imagery  from  the  collected 
data  following  the  experiment.  The  output  of  the  beamforming  algorithm  is  a  complex  value 
representing  the  magnitude  and  phase  of  the  backscattered  field  at  each  resolution  cell. 
The  remainder  of  the  signal  processing  follows  closely  techniques  used  in  conventional  wind 
profilers  in  which  the  Doppler  spectrum  of  the  echo  is  estimated  using  Fourier  transform 
or  pulse-pair  techniques.  Ground  clutter  and  any  spurious  signals  are  extracted  from  the 
spectra  prior  to  computing  the  first  three  moments:  total  power,  mean  velocity,  and  spectral 
width.  Volumetric  images  of  any  of  these  data  products  can  be  generated  as  often  as  ten 
times  per  minute. 

Figure  3  shows  a  subset  of  the  basic  TEP  data  products.  Figure  3a  shows  a  TEP 
image  of  the  backscatter  intensity  (C,^)  taken  at  Rock  Springs  in  August  of  1996.  The 
turbulence  data  shown  is  resolved  into  30-m  range  gates  between  the  altitudes  of  900 


(m/s)  dB  RELATIVE  TO  MEAN 


and  1290  meters.  Sequential  images  for  a  nine-minute  interval  appear  at  the  Web  site: 
http://acadia.ecs.umass.edu/html/tep-psu96.html.  Figure  3b  shows  the  radial  veloc¬ 
ity  sampled  concurrently  with  the  intensity  image  shown  in  figure  3a  Each  image  represents 
a  five-second  average  of  the  atmospheric  return. 

The  three-dimensional  fields  obtained  by  TEP  are  unique  compared  to  other  clear  air 
radars  because  the  beam  forming  process  allows  the  simultaneous  measurements  of  backscat- 
ter  and  velocity  at  every  pixel  [19].  Thus,  TEP  affords  a  synoptic  view  of  the  ABL  that 
is  not  currently  possible  with  other  remote  sensing  instruments.  However,  additional  re¬ 
search  is  needed  before  TEP  can  fulfill  its  promise  for  contributing  to  the  understanding 
of  ABL  turbulence.  Measurements  from  this  system  must  be  compared  to  and  velocity 
measurements  obtained  by  other  remote  sensing  systems  such  as  lidars  and  EM-CW  radars. 

During  the  last  year  of  the  URI  we  were  able  to  make  some  preliminary  comparisons  of 
TEP  data  with  PSU  LES  computations.  We  used  a  40  minute  record  of  TEP  data  from 
a  highly-convective,  afternoon  boundary  layer  (the  inversion  layer  height  Zi  =  1140  m,  the 
geostrophic  wind  Ug  =  0.9  m/s).  Our  PSU  colleagues  provided  us  with  and  velocity  data 
from  a  highly  convective  boundary  layer  {zi/L  =  —540,  Ug  =  1.0  m/s)  that  should  compare 
well  with  the  TEP  data. 

Qualitatively,  similar  features  can  be  seen  in  the  local  structure  of  the  TEP  and  LES 
data  sets.  Figure  4  shows  TEP  measured  horizontal  wind  vectors  overlayed  on  backscattered 
intensity  {CD  and  vertical  velocity  images,  respectively,  from  a  height  of  990  m  (z  =  O.STzj). 
There  is  a  convergence  of  the  horizontal  winds  on  to  the  area  of  high  relative  (7^,  corre¬ 
sponding  to  a  downdraft  feature  in  vertical  velocity.  Such  features  appear  frequently  in  the 
TEP  data  set,  and  always  correspond  to  a  small,  coherent  downdraft  feature  surrounded  by 
a  larger  updraft  feature.  LES  shows  similar  correspondance  between  converging  horizontal 
winds,  coherent  downdrafts  in  vertical  velocity,  and  locally  high  values.  Figure  5a  shows 
the  local  (7^  predictions  from  LES  sX  z  =  0.952:i  plotted  on  a  dB  scale  with  horizontal  veloc¬ 
ity  vectors  overlayed;  figure  5b  shows  the  corresponding  LES  vertical  velocity  image.  The 
qualitative  similarity  of  the  LES  predictions  and  the  TEP  measured  data  suggests  that  LES 
may  be  correctly  predicting  local  <7^  behavior. 

We  have  also  compared  TEP  and  LES  (7^  values  quantitatively.  By  taking  cross-wind 
slices  through  the  data  shown  in  Figure  3a  and  plotting  them  as  a  function  of  time,  larger 
volumes  are  created.  Such  volumes  assume  that  the  turbulence  is  frozen  in  space  and  time 
and  advected  through  the  radar  beam  by  the  mean  wind  [20].  Such  a  volume  is  shown  in 
figure  6,  and  provides  a  better  statistical  approximation  of  the  LES  domain. 

Peltier  and  Wyngaard  [16]  derived  statistics  of  local  structure-function  parameters  for 
humidity,  Cq,  humidity-temperature  correlation,  Ctq,  and  temperature,  C|.  for  a  moderately 
convective  boundary  layer  using  LES.  Through  the  relationship  in  [21], 

d  =  a^Cl  -h  2oiI5Ctq  + 

the  statistics  of  local  can  be  generated  over  the  LES  domain.  A  “variability  index”,  or 
normalized  variance,  can  be  defined, 

p  variance  of  Cl 
"  squared  mean  of  C^ 

which  is  computed  over  horizontal  slices  throughout  the  LES  volume.  Using  equation  1,  we 
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Figure  4:  TEP  convergence  of  horizontal  winds  on  an  area  of  high  local  C^- 
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Figure  5:  LES  convergence  of  horizontal  winds  on  an  area  of  high  local 
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have  compared  TEP  and  LES  statistics.  Figure  7  shows  that  comparison,  with  each  data  set 
scaled  to  the  boundary  layer  height  Zi  .  The  two  curves  exhibit  a  similar  trend  from  0.55  Zi 
upward,  although  it  appears  as  though  there  is  a  small  offset  between  the  two  curves.  One 
hypothesis  we  have  for  that  offset  is  that  the  beam  efficiency  of  the  TEP  system  is  lower 
than  we  have  anticipated.  Future  efforts  to  improve  the  reliability  and  calibration  of  TEP 
should  improve  the  beam  efficiency  and  test  that  hypothesis.  Other  sources  of  error  could 
include  the  effects  of  the  field  site  terrain  and  the  length  of  the  data  set,  and  the  experiments 
as  part  of  this  program  will  study  those  effects. 

Below  0.55  Xj,  the  for  the  TEP  data  is  quite  large.  We  believe  that  those  effects  are 
due  to  biological  targets,  such  as  bugs  and  birds,  in  the  lower  boundary  layer.  Such  targets 
are  very  bright  in  radar  backscatter  and  are  quite  intermittant,  producing  a  large  F„.  The 
deployment  of  a  millimeter  wave  radar  along  with  the  TEP  system  would  allow  MIRSL 
researchers  to  identify  and  remove  from  the  data  non-atmospheric  targets. 
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